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ABSTRACT 

We present results of -'^^CO (1~0) and ^'^CO (1-0) observations of the northeastern spiral arm 
segment of IC 342 with a ~ 50pc resolution carried out with the Nobeyama Millimeter Array. Zero- 
spacing components were recovered by combining with the existing data taken with the Nobeyama 45m 
telescope. The objective of this study is to investigate the variation of cloud properties across the spiral 
arm with a resolution comparable to the size of giant molecular clouds (GMCs). The observations 
cover a 1 kpc x 1.5 kpc region located ~ 2 kpc away from the galactic center, where a giant molecular 
association is located at trailing side and associated star forming regions at leading side. The spiral 
arm segment was resolved into a number of clouds whose size, temperature and surface mass density 
are comparable to typical GMCs in the Galaxy. Twenty-six clouds were identified from the combined 
data cube and the identified clouds followed the line width-size relation of the Galactic GMCs. The 
identified GMCs were divided into two categories according to whether they are associated with star 
formation activity or not. Comparison between both categories indicated that the active GMCs are 
more massive, have smaller line width, and are closer to virial equilibrium compared to the quiescent 
GMCs. These variations of the GMC properties suggest that dissipation of excess kinetic energy of 
CMC is a required condition for the onset of massive star formation. 

Subject headings: galaxies: individual (IC342) - galaxies: ISM - galaxies: spiral - ISM:clouds- 
ISM: molecules 



1. INTRODUCTION 

The formation of massive star is one of the fundamen- 
tal processes driving secular evolution of spiral galaxies. 
The rate of massive star formation in galaxies is consid- 
ered to be regulated by the formation and the evolution 
of giant molecular clouds (GMCs), which are progeni- 
tors of massive stars. However, the exact processes which 
drive the formation and the evolution of GMC to initiate 
massive star formation is still unclear. 
GMCs are large cloud complexes with mass of 10^ to 
several times of 10 ^ Mp) and size of 20 - lOOpc (e.g., 
Sanders et all 119851: iSolomon etall 119871: iDame et alj 



Strong scaling relations between the size and line width 
of GMCs wer e found both in the Galactic clouds (e.g. , 
LarsonI llQSlt IDame et al.l 119861 : ISolomon etall 119871 : 



1986f ). GMCs comprise roughly 80 % of the total niolec- 
ular mass in the i nner Galaxy ([Sanders et al.l 119851 : 
ISolomon etHlllQSTt ). Ahhough the mean volume den- 
sity of GMCs is as low as ~ 100 cm~^, there is strong 
density contrasts inside the GMCs, which make local vol- 
ume density higher than 10^ cm~^, which is thought as 
critical density f or the onset of the star formation (e.g., 
lElmegreenI I2002D . As formation of dense clump is es- 
sential in forming stars, interpretation of the physical 
process which make a density contrast inside GMC is 
crucial. 
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Hever fc BruntI I2004D and in the extragalactic clouds 
(e.g.. iRosolowskv et al.l[200l iBolatto et al.ll2003D . The 

scaling relation holds not only across clouds but also 
within clouds, suggesting that entire structures of GMCs 
are strongly governed by su personic turbulent motion 
(e.g. lLarsonl[l98lt iHever fcBru nt 20041 . The supersonic 
turbulence act in two ways: while it supports the GMCs 
against global gravitational collapse, it also plays an im- 
portant role in forming density stru ctures within each 
CMC (turbulent fragmentation, e.g., iVazguez-Semadenil 
119941: [Padoan fe Nordlund 2002). There is an argu- 
ment that, as a consequence of the fragmentation pro- 
cess, dense cloud cores, which finally collapse into 
stars, are formed inside GMCs (e.g., Elmegrecn 200% 
iMac Low fc Kle ssen^200"l iKrumholz fc McKeei2005l , al- 
thou g:h additional process, such as competitiv e accre- 
tion (jBonnell et al.l l200ll: iBonnell fc Bati [20061 ) . might 
be also required in massive star formation. An investiga- 
tion of the dynamical properties of GMCs is important 
in discussing the formation and evolution processes of 
GMCs. 

Recent molecular cloud studies in the Galaxy indi- 
cated that GMCs which lack associated star forma- 
tion activity are not as unusual as once thought (e.g., 
IChiar et al.lir99l . A comparison of cloud properties be- 
tween star-forming GMCs and quiescent GMCs is impor- 
tant in addressing what re ally control s the star forma- 
tion activity in GMCs. Williams et al.l ([19941 compared 
cloud properties between a qu iescent cloud (G216-2.6, 
known as Maddalena's cloud, iMaddalena fc ThaddeusI 
1985f) and an active clou d (the Rosette molecular cloud. 
Blitz fc ThaddeusI [l980l ). They showed that the line 



width is larger in the Maddalena's cloud than in the 
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Rosette cloud over every size scale. However, distance 
degeneracy in observations of Galactic clouds prevents 
further investigation of the relation between cloud prop- 
erties and the star formation activity. To overcome this 
limitation, observations of molecular clouds in external 
galaxies are crucial. 

Early millimeter molecular-line observations in the exter- 
nal grand-design spiral galaxies often found large associ- 
ations of molecular clo uds (e.g.. lVogel et al.lll988l : iRandl 
ll995l : lRand et al.lll999| ). Mass of those clouds were larger 
than GMCs by over an order of magnitude (greater than 
several times of 10''' M(^) and termed as giant molecular 
association (GMA) (jVogel et all 119881 ). Distributions of 
those GMAs were found to be coherently offset from as- 
sociated Hn regions and interpreted as representative of 
the time delay between the formation of those massive 
clouds and t he onset of the subsequent massive star for- 
mation fe.g.. lRandlll993D . By resolving the spatial offsets 
down to the scale of GMCs, it should be able to investi- 
gate the evolution sequence of GMCs which leads to the 
massive star formation. 

We report the results of the observations of the nearby 
galaxy IC 342 carried out with the Nobeyama Millimeter 
Array (NMA) in ^^CO (1-0) and ^^CO (1-0) lines. The 
observed field covers a 1 kpc x 1.5 kpc region located in 
the northeastern spiral arm segments, 2kpc away from 
the galactic center. At the site, spatial offsets between 
molecular clouds and star forming regions were previ- 
ously MentifedbytheNobey 45m telescope observa- 
tion (|Hirota et al.ll201(il) . The prime target of this study 
is to investigate variation of cloud properties across a spi- 
ral arm in IC 342 with a spa tial resolution comparable 
to th e size of GMCs (- 50pc: iSanders et all figSS: Bliti 
[19931) . 

IC 342 is a nearby spi ral galaxy classified as SAB(rs)cd 
()de Vaucouleurs et al.l lT991). The central region of IC 
342 harbors starburst activity fe.g.. iBecklin et aLlll980( ) 
and because of the strong millimeter and submillime- 
ter emission from the region, the central region (within 
~ 500pc) of the galaxy has been studied in detail. 
The distributions of molecular c l ouds is indicative of 
bar p otential (e.g., Ilshizuki et al.l [l990l : [Turner &: Hurti 
Il992t ). and molecular g ases are likely to be excited by 
starburst activity ( e.g., 'Wall & Jaff^'1990'; 'Eckar t et al. 
1990; Turn er et alT l993: Schulz ct al. 2001; Meie r et al. 
2000t Israel & Baas 2003). Studies of molecular chem- 
istry indicated the chemistry of the cl ouds in the cen- 
tral region are also sub ject to shocks ([Meier fc Turner! 
I2005t lUsero et~an I2006D presumably driven by the bar 
kinematics. Supply of material driven by the bar kine- 
matics is likely to be regulated by negative mec hanical 
feedback from the starburst activity (Schinncre r et ahl 
I2008D . Beyond the radius of 500pc, the distribution 
of the molecular gas disk of IC 342 was studied with 
the single dish telescope observation mostly in ^^CO (1- 
0) f Sage fc Solomon"1991l: iCrosthwaite et al.ll2001l : iSato I 
[2006; Kuno ct al. 200 7,). In part i cular^ the ^^CO (1- 
0) map presented in iKuno et al.l ([2007| ) reveal promi- 
nent galactic structures such as a bar and spiral arms. 
These structures ar e also visible in near-infrared images 
([Jarrett et al.ll2003[ ). suggesting the existence of moder- 
ately stron g density waves i n the galaxy. The proximity 
(3.3 Mpc, iSaha et all I2002D of the galaxy provides an 



opportunity to investigate the influence of density waves 
exerted on molecular clouds. 

2. OBSERVATION AND DATA REDUCTION 

2.1. 12(70 (1-0) and CO (1-0) observation 

Aperture synthesis observations of the northeastern 
spiral arm segments of IC 342 in 1^00(1-0) and ^^CO 
(1-0) were carried out with the NMA. The NMA con- 
sists of six 10m antennas which provide fleld-of-view size 
of - 59" and - 62" at the rest frequency of ^^CO (1-0) 
and i^CO (1-0), respectively. Three antenna conflgura- 
tions (AB-f C-l-D) were used to sample visibility in the 
uv range of '-^ 4 fcA to ~ 130 k\. The observed region 
covers a GMA located on the northeastern spiral arm seg- 
ments, where spatial offsets among ^^ CO (1-0), ^'^CO (1 - 
0) and star formation tracers exist ([Hirota et al.l 120 iol ). 
To cover the region of interest, two pointings separated 
from each other by 30" were set (Figured]). The ob- 
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Figure 1. Velocity-integrated image of ^^CO (1-0) data taken 
from Sato (2006 ) (contour) compared with 24/jm imag e taken with 
the M ultiband Imaging Photometer for Spitzer (MIPS: fRieke et al.l 
l200i) camera. Contour levels are 10, 15, 20, 25, 30, 35, 40, 50, 
70, 100, 140, and 190 K km s"!. Pointing centers of the NMA 
observations are indicated with crosses. Dashed circles with size of 
60" indicate approximate size of the NMA field of view (~ 59" for 
i^CO (1-0) and ~ 62" for "CO (1-0)). The peak position of the 
GMA seen in the ^^CO image is indicated with a diamond. 



servations were carried out from 2005 November to 2007 
April for ^^CO (1-0) and from 2006 November to 2007 
March for ^■^CO (1-0). System noise temperatures (in 
single side band) were 500-lOOOK for ^^CQ (1-0) obser- 
vations and 400-700K for "C O (1-0) observations The 
Ultra- Wide-Band-Correlator ([Okumura et al.|[2OOO0 con- 
figured to cover 256MHz bandwidth with 256 channels 
was used as the backend. Window function for the spec- 
trometer was set to Banning smooth function and the 
resultant frequency resolution was 2MHz. Two point- 
ings were observed alternately every 8 minutes to attain 
a uniform Mw-coverage. B0355-I-508 was observed once in 
every ~ 20 minutes as a gain calibrator and 3C273 was 
observed to determine the pass-band. 
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Acquired raw visibility data we re calibrated us i ng th e 
software package UVPROC-II (|Tsutsumi et alj I1997D . 
Though the NMA data were combined with the data ob- 
served with the 45ni telescope afterward, map image of 
the CO lines were first made without combining single- 
dish data to check the basic quality of the data. Imaging 
and deconvolution were made following the standard pro- 
cedures imp lemented in the software package MIRIAD 
(jSault et al.| [T995). Resolution of the resultant ^^CO 
data cube was 3".l x 2". 5 in spatial directions and 5.2 
km s~^ in velocity direction. Typical rms noise was ~ 24 
mJy beam~^ within each channel. For ^^CO (1-0) data 
cube, spatial and velocity resolutions were 3". 88 x 3". 53 
and 

~ 5.4km s~^, respectively, and typical rms noise was 
~ 11 mJy within each channel. The total amount of 
the flux detected with the NMA is 8.7 x 10^ Jy km 
s-i for I^CO (1-0) and ~ 4.0 Jy km s^^ for ^^CO (1- 
0). While, on the other hand, the total flux emitted 
from the observed region was estimated from the 45m 
data and was ~ 6.1 x 10^ Jy km s'^ for ^^CO (1-0) and 
~ 8.7 X 10^ Jy km s'^ for ^^CO (1-0). The large amount 
of the 'missing' flux discrepancies indicates the presence 
of extended structures whose spatial frequency is lower 
than ~ 4.5 fcA, which is the minimum spatial frequency 
available with the NMA. 

2.2. NMA and 45m Combined Data: Short Spacing 
Correction 

To recover the missing extended flux, the NMA data 
were combined with the uv data generated from the 45m 
data to fill in the central 'h ole' in the Mz ;-coverage. ^^CO 
(1-0) data were taken fr o mlSatol poM ) and "CO (1-0) 
data from iHirota et al.l (|2010l ). Combin ing procedures 
foll ow the method descr ibed in Takaku wa et al.l (|2003f ) 
and iKurono et al.l ()2009l ) . We will denote outlines of the 
procedures hereafter. The 45m data were deconvolved 
with a Gaussian pattern using a Wiener filter and mul- 
tiplied with the NMA primary beam pattern to simulate 
the NMA observation. Visibility data were generated 
from the 45ni data by performing Fourier transforma- 
tion on the deconvolved 45m data. On generating vis- 
ibility data, (m, u)-sampling points were uniformly dis- 
tributed within the radius of 4.5 kA, which corresponds 
to the 'central hole' of the NMA {u, u)-coverage. Both 
the NMA and the 45m visibility data sets were merged to 
create a single visibility data set and imaged with the in- 
verse Fourier transformation. " Clean" deconvolution was 
performed with the MOSSDI task implemented in the 
MIRIAD. Finally, the images were corrected for primary- 
beam attenuation by dividing with the gain distribution 
estimated with the standard MIRIAD task (MOSSEN). 
On performing the inverse Fourier transformation, 
weights attached to the 45m and the NMA visibility de- 
termine the shape of the synthesized beam. If larger 
weights are attached to the 45m data, sensitivity to the 
low spatial frequency components in increased. How- 
ever, if the 45m data are overweighted, the shape of the 
synthesized beam gets close to the 45m beam and thus 
resolution gets worse. 

In practice, the relative weights between the 45m data 
and the NMA data are determined by the number of the 
(u, w)-sampling points generated from the 45m data and 
the integration time attached to each visibility sample 



for the 45m data. An adequate choice of the parameters 
should attain the resolution comparable to the original 
NM A data and th e complete recovery of the missing- 
flux (jKurono et al.l 20091. We generated several sets of 
visibility data for the 45m data with different integration 
time attached to each visibility data sample. The number 
of visibility data points were fixed since only the product 
of the both parameters is important. Best weight param- 
eter, which satisfies the requirement of the least amount 
of the missing flux and the least broadening of the beam, 
was selected from the trial data sets. 
Properties of the flnal combined data cubes are listed in 
Table [1] Spatial resolutions of the combined data cubes 
are almost identical to the original data cube. The total 
flux within the observed fields are also almost identical 
to the single-dish data within the ~ 5 % level, which is 
well below calibration errors. 



3. RESULTS 

3.1. Molecular gas distribution: ^^CO (1-0) 

The observed region contains a GM A which was visi- 
ble with the previo us CO observations ()Kuno et al.ll2007l : 
IHirota et al.ll2010[ ). The NMA-h45m combined data of- 
fer an opportunity to investigate the internal-structure 
of the GMA with a spatial resolution comparable to the 
size of GMCs. Figure [2] shows the velocity-integrated in- 
tensity image of the combined ^^CO(l-O) data. Several 
discrete sources with sizes of 50-lOOpc were found. For 
convenience, three brightest ^^CO sources were marked 
with crosses in the figure (Points 1-3). Indices were al- 
located in order of the ^^CO intensity. Typical bright- 
ness temperatures (Tb) at the peaks are 2-4 K above the 
cosmic background temperature, with a maximum tem- 
perature of 5.2 K at Point 1. Although these values are 
lower than the peak CO temperatures found in Galactic 
GMCs (~ 20K), it is comparable with the mean temper- 
ature averaged over f ull extent of cloud. For example, 
iSakamoto et al.l ()1994D derived the spectrum of ^^CO (1- 
0) and ^^CO (2-1) averaged over almost entire extent of 
the Orion A and B clouds and presented ^^CO (2-1) peak 
intensities of 1.5K and 2. IK, and CO (2-1) / CO (1-0) 
ratio of 0.75 and 0.62. These values correspond to ^^CO 
(1-0) temperature of about 2. IK and 3.4K for both the 
clouds. Comparable brightness temperatures found with 
the NMA-f 45m data suggests that clouds seen with our 
observations may have similar nature with the Galactic 
clouds. 

The characteristics of ^^CO distribution on both the 
southern and the northern sides of the observed fields 
differ apparently: while it is concentrated in the narrow 
ridge structure on the southern side, it is more extended 
on the northern side. On the southern side, the molecular 
ridge consists of several discrete sources including Points 
1 and 2. In a contrasting situation, on the northern side, 
the cloud distributions are more extended compared to 
the southern side. The GMA seen as a single peak in 
the previous 45m image was resolved into several minor 
peaks around Point 3. 

The amount of molecular gas mass is estimated from the 
CO integrated intensity by applying a "standard" CO- 
H2 conversion factor of Xqq = 2 x 10^° cm~^ (K km 
s-i)^i (|Strong fc MattoxlligM iDame et al.llMl . The 
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Table 1 

Parameters of the Combined CO Data Cubes 



Line 


Beam Size 


Beam Position Angle 


Velocity resolution 


rms Noise 




(arcsec) 


(deg) 


(km s^^) 


(mjy beam^^) 


12C0 (1-0) 


3". 20 X 2". 58 


-72° 


5.2 


23 


13C0 (1-0) 


3". 85 X 3". 53 


-69° 


5.4 


10 



a) '^CO (NMA+45m) 



IJJy/beam km/^ 



b) 'XO (NMA+45m) 
0.0 0.5 



^ I JJy/beam km/s] 



68 40' 
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Figure 2. ^^CO (1-0) and ^'^CO (1-0) integrated intensity images of the northeastern spiral arm segments of IC 342. Loeations of the 
three brightest ^^CO sources are indicated with crosses (Points 1-3). (a) ^^CO (1-0) integrated intensity map of the NMA+ 45m combined 
data. Contour levels are 2, 4, 6, 8, 12, and 20 times 0.45 Jy beam^^km s~^. Dashed circles indicate the approximate size of the NMA field 
of view (^ 60"). The attenuation due to the primary beam pattern is corrected, (b) Same as (a), but for the ^^CO (1-0) data. Contour 
levels are 1, 3, 5, 7, and 9 times 0.15 Jy beam~^km s~^. 



surface mass density of H2 is calculated by 



3.21xcosixf — (1) 



where i is the inc l inatio n angle of the galaxy (i ~ 31°; 
iCrosthwaite et all l2000( ) . The surface mass density of 
molecular gas corrected for He and other heavy elements 
is derived by 

I],„o1-1.36xSh,. (2) 

Highest peak surface molecular gas mass density within 
the observed region is ~ 23OM0 pc~^ found at Point 
1. Typical surface mass densities in other local 

peaks are 100 2OOM0 pc^^. These are consistent 

with the value 100M(T) pc~ ^) found in the Galac- 



tic clouds (e.g., iSolomon et al. _1987: Blitz 1993 il and 
several nearby galaxies (e.g.. iRosolowskv et al.l 120031 : 
[Rosolowskv fc Blitz 20051) . The size, the temperature 
and the surface mass density of the discrete sources found 
in the ^^CO map imply that the NMA-|-45m observations 
successfully resolved the molecular spiral arm into indi- 
vidual GMCs. 

3.2. Molecular gas distribution: ^^CO (1-0) 



Figure [2][b) shows an integrated intensity map of the 
combined ^"^CO (1-0) data. To improve the signal-to- 
noise ratio (S/N) of the integrated image, ^^CO data 
were masked before calculating the integrated intensity. 
The mask data were made according to the following pro- 
cedures. First, the ^^CO data cube was smoothed and 
re-sampled to share the same resolution and pixel coor- 
dinates with the ^'^CO data cube. Mask data was made 
from the smoothed ^^CO dat a cube according to the fo l- 
lowing two-step procedures f Rosolo wskv &: Blitz! I2005D : 
first, pixels above 4cr are adopted as 'kernel' mask and 
next, all the pixels above 2cr and connected with the 'ker- 
nel' are included within the mask. Isolated masks smaller 
than the beam size were discarded. The applicapability 
of this mask to the "'^^CO data relies on the following is- 
sue: as the critical volume density for the excitation of 
^^CO line 10^ cm~^) is much lower than that of ^^CO 
line (~ 10'^ cm~^), the extent of the ^^CO emitting vol- 
umes should be enclosed within that of ^^CO emission 
regions. To check the validity of the mask, intensity his- 
togram of ^^CO for the residual pixels outside the mask 
was made. The shape of the residual ^^CO intensity his- 
togram was symmetric about zero and the central part 
of the histogram could be fitted with the Gaussian dis- 
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tribution with a sigma of ~ 11 mJybeam" , which is 
consistent with the rms noise of the original ^^CO data, 
suggesting that almost all the ^'^CO emission were en- 
closed within the mask. Finally, the mask was applied 
to the original ^'^CO data cube and integrated intensity 
was calculated. 

While ^'^CO peaks on the southern side are aligned 
with ^^CO peaks, ^"^CO peaks on the northern sides are 
largely deviated from ^^CO peaks. In particular, the 
^^CO sources located around Point 3 (GMA center) lack 
clear counterpart in ^■^CO image. As in the case of ^^CO, 
the ^'^CO distribution on the northern side is spatially 
extended compared to the southern side. 

3.3. Comparison with Multi-wavelength Data 
3.3.1. Star formation tracers 

The ^^CO image was compared with Ha, 8/xm, and 
24/xm images to see the spatial relation between the 
distributions of molecular clouds and star forming re- 
gions. Monochromatic Ha image data were provided by 
Hernandez et al.l(|2005D . Both the 8/im image taken with 
the Infrared Array Camera (IRAC; iFazio et all 120041) 
and the 24/im image taken with the Multiband Imag- 
ing Photometer ([Rickc ct al. 2003) were retrieved from 
the Spitzer archive. Resolution of the Spitzer images is 
~ 1".2 for the 8/Ltm image and ~ 6" for the 24/Ltm image. 
Figure |31[a)-(c) show the Ha image, the 8/Ltm image, and 
the 24/im image respectively. All the images are overlaid 
with the i^CO (1-0) image. 

It is apparent that while the clouds on the southern side 
are closely associated with their neighboring star form- 
ing regions most of the clouds on the northern side seem 
to lack associated star forming regions. In particular, 
around Point 3 (center of the GMA), little star forma- 
tion activities are seen seen both in the Ha and the mid- 
infrared images. As we have seen before, the molecular 
gas distribution around the center of the GMA is much 
more smooth compared to in the southern side where 
molecular clouds are concentrated to form the narrow 
ridge. To see whether such difference of the distribution 
and star formation activity in both sides of the observed 
region is related to the properties of clouds, cloud prop- 
erties will be examined in the latter (Section 14. 1|) . 
Previous CO observations made with coarse spatial res- 
olutions (300-1000 pc) often found the well-ordered spa- 
tial offsets between the molecular spiral ar m and star- 
forming re gions in other spiral galaxies fe.g.. lVogel et al.l 
[1988 : Rand et al.l 11999.) . Also in IC 342, with the spa- 
tial resolution of ~ 320p c, the GMA is seen offset from 
the star-forming regions (|Hirota et al.l IMcil) . While on 
the other hand, seen with the 50pc resolution here, the 
separations between the clouds and the associated star- 
forming regions are small (mostly below the beam size). 
Recent observations of the gra nd design spiral g alaxy 
M51 also indicate similar result (lEgusa et al.ll201lf ). It is 
suggested that large GMA found with the previous coarse 
resolutions is a mixture of the both kinds of clouds, which 
are clouds associated with and not with star forming re- 
gion, as in IC 342. 

3.3.2. Hi 

Figure |3]d) shows the comparison between the ^^CO 
image and the H I image. H I data were retrieved from 



the Very Large Array archive and reduced with the As- 
tronomical Image Processing System. The resultant H I 
data cube had a spatial resolution of 23". 2 x 20". 1 and 
typical rms noise of ~ 0.7 mJy beam~^ for each channel. 
The surface mass density of H I is calculated under the 
assumption of optically thin H I emission. Typical sur- 
face density of H I over the observed field is ^ 2 4 

Mq pc~^ and is smaller than that of molecular gas by 
an order of magnitude. Though the spatial resolution of 
the H I image is ~ 20" and worse than the NMA image, 
it is apparent that there exists a hole at the center of the 
GMA. Around the H I hole, distribution of H I roughly 
coincides with that of H II regions. The coincidence of 
H I with H II regions indicates that the H I clouds around 
the GMA might produced by the dissociation of H2 by 
the UV radiation. 

3.4. Cloud identification and basic cloud properties 

To investigate the properties of molecular clouds across 
the spiral arm, decomposition of individual cloud emis- 
sion from the combined ^^CO (1-0) data was attem pted. 
The CLUMPFIND algorithm (|Williams et al.llT 99?) with 
some fine-tuning of p aram eters, which were proposed by 
iRosolowskv fc Blitz! ([2005I ). was utilized. Original code 
of the CLUMPFIND algorithm was targeted at identify- 
ing Galactic molecular cloud data observed with single- 
dish telescopes, in which beam size of the observation is 
comparable to the pixel size and S/N greatly differs from 
typical extragalactic observations. Rosolowskv fc Blitl 
((2005j) performed cloud decomposition from the interfer- 
ometric data of M64, which is a molecular rich galaxy 
with a similar distance to IC 342. They proposed some 
modifications to be made to the original CLUMPFIND 
algorithm in applying the algorithm to extragalactic ob- 
servation data. As the spati al resolution and the S/N o f 
our data resembles that of IRosolowskv fc Blitz! (|2005D . 
we adopted some of the modifications proposed by the 
author. 

The modifications proposed bv IRosolowskv fc Blitd 

(j2005f ) consists of two points: one is to alter the man- 
ner of locating significant cloud peaks and another is to 
change distance metric along the velocity axis to account 
for the oversampling in spatial direction. Among these 
two modifications, we had adopted only the latter one. 
In the following, we will briefiy mention about the mod- 
ification we had adopted and not adopted. First, we had 
adopted the change of distance metric. As the spatial res- 
olution of the combined data is ~ 3" and the cell size of 
the data is 0.4", the data are heavily oversampled in the 
spatial direction compared to the velocity direction, in 
which the resolution is ^ 5.2 km s~^ and the pixel size is 
^ 2.6 km s~^. To account for such oversampling, the dis- 
tance metric used in the parti tioning cloud boundary was 
modified (see Appendix A of IRosolowskv fc Blitz 2005) . 
Second, we tested the modified procedure of finding sig- 
nificant local maxima. While the original CLUMPFIND 
code searches and partitions data cube into each clump 
by a fixed intensity interval (usually taken as 2(t), the 
method proposed by.Rosolowsky fc Blitz (2005) searches 
and extends clumps by every 0.5cr but adopt 3cr thresh- 
old in discarding 'false' local maxima. We tested the 
both algorithms and found that the results of the cloud 
participating are similar to each other, except for clouds 
around Point 3 (the GMA center). While the original 
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Figure 3. Contour map of the velocity-integrated ^^CO (1-0) intensity image overlaid on the pseudocolor scale representations of (a) 
the H« image, (b) the S/im image, (c) the 24/im image, and (c) the H I image, respectively. Contour levels are 1, 3, 5, 7, and 9 times 1.0 
Jy beam~^ km s~^. Dashed circles indicate the approximate size of the NMA field of view (^ 60"). 



CLUMPFIND code produced three clouds around the 
point, the modified CLUMPFIND algorithm bundled the 
three clouds into a single large cloud complex. However, 
close examination of the results indicated that the three 
clouds bunched by the modified CLUMPFIND algorithm 
seem to be individual entity because of the different cen- 
ter velocity. Because of this, we did not adopt the mod- 
ified procedure for finding local maxima in participating 
the data. 

Identification of clouds from the combined NMA-f 45m 
data was performed with the following steps. A data 
cube containing S/N value as a pixel value was gen- 
erated by dividing the data cube before primary beam 
correction with rms noise level. The CLUMPFIND al- 



gorithm ([Williams et al.lll994D with minor modifications 
mentioned above was applied to the S/N data cube with 
a 2a increment level and 2cr lowest cutoff level. After par- 
titioning the data cube into each clump, clumps located 
completely outside of the primary beam and clumps with 
peak temperature below 7a were excluded from the anal- 
ysis. The 7a cutoff was set by the fact that over this level, 
both the algorithms tried here produced similar results, 
as have noted in the previous. Finally, 26 GMCs were 
identified from the ^^CO data cube. Figure [4] shows lo- 
cations and approximate extent of the identified clumps 
projected on spatial directions. Figure [5] shows channel 
maps of the selected regions from the data cube with 
peak positions of the identified clumps overlaid, to show 
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Figure 4. Locations of the identified GMCs overlaid on the ve- 
locity integrated ^'^CO (1-0) image. Contour levels are the same 
as in Figure [2] Positions of the identified clouds are indicated 
with crosses. Ellipses indicate the FWHM sizes of the intensity 
distribution for each cloud. The FWHM sizes and position an- 
gle of the ellipses were determined by calculating the eigenvec- 
tors of the i ntensity- weighted covariancc matrix for each cloud (see 
et al.l f2006; Rosolowsky Sz Lerov. .2006> ). Clouds IDs are also 
indicated on the right of the crosses. 



how the identification of the clumps was executed. 
Properties of the identified clouds, namely, mean posi- 
tion, mean velocity, cloud radius (R), FWHM velocity 
width (Ay), and CO luminosity (L(^^CO)) were mea- 
sured with the data cube corrected for primary beam 
attenuation. The mean position and the mean velocity 
of the clouds were taken as the first moment of the emis- 
sion. The effective radius of the identified clouds was 
calculated by 



2 ^l/4 

cam,inaj 



ocam.inin/ 



(3) 



where 3A/^/tt is a factor for conve rting the rms size int o 
the radius of a spherical cloud (jSolomon et al.l [T987D . 
(Tx, (Jy is the second moment of the intensity distribu- 



tion in spatial directions, and (Tb 



cam,maj 



and (Jbc 



is 



the rms size of the observed beam in major and minor 
axis directions, respectively. As the rms sizes {ax^ (Jy) 
are underestimated because of the 2a clipping level for 
each cloud boundary, correction should be made. Of- 
ten used method in this case is a Gaussian correction 
which assumes a Gaussian profile of the cloud emis- 
sion and boosts the measured properties with the fac- 
tor determined by the ratio of the peak temp erature to 
the truncation level ( Oka et al. 2001; Bolatto et al.ll2003l: 
iRosolowskv fc Blitzl i2005t iRosolowskvi ,2007.) . The cor- 
re ction was made following th e analytic expression given 
bv lRosolowskv fc Blitd ()2005D . which boosts the rms size 
by factor of at most 1.2 for our data. Line width of 
the clouds was taken as the FWHM size of the cloud 
profile, also corrected with the Gaussian correction, and 
corrected for resolution bias by subtracting the velocity 
resolution 2.6 km s~^) in quadrature. CO luminos- 
ity was taken as the summation of the emission within 



each cloud boundary corrected with the Gaussian correc- 
tion. Luminosity based cloud mass (Mco) was derived 
from the CO luminosity by applying the standard con- 
version factor noted in the previous section fSection l3.1l) . 
The virial mass of a spherical cl oud with density profile 
p oc r~" is written as (e.g., MacLaren et al. li 



= 126- 



2n 



3 — ?i 



Ay 

km s~ 



(4) 



The density profile of n = 1 was assumed in calculating 
the virial mass. 

Large scale kinematics of the galactic disk such as the 
galactic rotation and streaming motions may bias the 
measured line width and hence the virial mass. To check 
the possible infiuence of this, the velocity shear across 
each cloud was estimated by 



AKLar 



(5) 



where Xi, yi is the position of each pixel in the cloud, 
Viot{xi^ yi) is the circular rotational velocity at point [xi, 
yi), and Wrot(a;o, J/o) is the circular rotational velocity at 
the mean position of the cloud (xq, yo)- An estimated 
amount of AT4hcar was smaller than the measured line 
widths by nearly an order of magnitude and turns out to 
be negligible for our measurement. 

Table [2] denotes the basic properties of the identified 
molecular clouds. The summation of the mass of the 
identified clouds is ~ 4.4 xIO^A/q, which comprises ~ 
60% of the total molecular mass within the observed re- 
gion. The relations between the basic properties of the 
identified clouds are indicated in Figure [6Ka)-(d). As the 
Gaussian correction is the largest source of uncertainty 
for the measurement of cloud properties, uncorrected val- 
ues are also indicated for comparison. 

3.5. Spatial Relation Between Star Formation Activity 
and Spiral Arm 

The comparison between the ^^CO image and the star 
formation tracers indicated that while some of the clouds 
seem to lack associated star forming activity, the rest of 
the clouds are closely associated with accompanying star 
forming regions. To check whether the properties of the 
clouds change with associated star forming activity, the 
identified GMCs were divided into two groups according 
to their star formation activity. 

The Ha and the mid-infrared images were examined 
for the classification. Although the 24/L(m band im- 
age is preferable in tracing the star formation rate 
(Calzetti et al. 2005, 2007), coarse resolution of the 24/im 
image {'^ 6") hampers identification of the correspon- 
dence between star-forming regions and the identified 
clouds. Because of this, the star- forming regions were 
identified from the IRAC 8/im image because of its higher 
resolution (~ 2") compared to the 24/im image. There is 
two possible bias in tracing star formation with the 8/im 
image: although the 8//m band is dominated by poly- 
cyclic aromatic hydrocarbon (PAH) emission mainly ex- 
cited by photon-dominated region around star-forming 
regions, there is little but certain amount of stellar con- 
tribution to the 8/im band (e.g.. iHelou et al]|2004D and 
the scaling between the star formation rate and PAH 
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RA (J2000) 

Figure 5. Channel maps of ^^CO(l-O) line emission for the selected spatial and velocity range of the data. Peak positions and IDs of the 
identified clouds are indicated with crosses and associated labels. Contour levels are 4, 6, 8, 10, 12, 14, 18, and 22 times 23 mjy beam^^, 
respectively. Label on the top left corner of each panel denotes the corresponding center velocity (Vlsr in km s"'^) for each channel. 



emission is slightly deviated from linear possibly because 
of the presence of diffuse PAH emission excit ed by ambi- 
ent radiation field (e.g.. ICalzetti et aLll2005D . However, 
as we are interested in locating star-forming regions and 
not in measuring the exact star formation rate, this point 
is not a severe deficit. 

As intensity distribution of star-forming regions seems 
to have distinct outline, t wo-dimensional versi on of the 
CLUMPFIND algorithm (jWilliams et aLllTOOl was uti- 
lized in defying the boundary of each star-forming region. 
Partitioning of the image was done by every 1.4 MJy 
str-i (~ 2a) step down to 7 MJy str'^ level 10a). 
The lowest boundary level was determined to keep away 
from diffuse 8/im emission which shall not be related with 
current massive star formation activity. The distribution 
of the identified star forming regions was also confirmed 
by comparing with the Ha image. 

The GMCs were divided into two group according to 
whether they are associated with the identified star form- 
ing regions or not. If a GMC overlaps with the identified 
star forming region and separation between the center of 
the GMC and that of the associating star forming region 
is within the beam size, the GMC is categorized as being 
"associated with star formation". Those "GMCs with 
H II region" are hereafter referred to as wHH" GMCs. 
On the other hand, the rest of the GMCs were cate- 
gorized as "GMCs without H II region" and termed as 



"woHH" GMCs. For most of the "wHH" GMCs, associ- 
ated star forming regions are seen in both the 8/zm and 
the Ha images, except for a GMC located near Point 
1. At Point 1, there exists an accompanying bright 8/im 
source but there is no clear counterpart in the Ha image. 
To qualitatively ensure the classification, 8/im fiux was 
measured within the boundary of each cloud. Figure [7] 
shows a histogram of the measured 8/xm flux. Although 
there is a "wHH" GMC which shows exceptionally low 
8/im fiux compared to other "wHH" GMCs (cloud-26), 
most of the "wHH" GMCs show higher 8/im fiux com- 
pared to the "woHH" GMCs. The low 8/im flux for 
the cloud-26 is due to the fact that the cloud is offset 
from the associated H II region and the Sfiui flux is mea- 
sured only within the cloud boundary. For the rest of 
the "wHH" clouds, offsets between each cloud and the 
associated star-forming region are not so large compared 
to the cloud-26. The 8/im flux histogram seems to ensure 
the classification made here. 

Figure |8] shows the spatial distribution of the classified 
GMCs. As in other grand design spiral galaxies, the 
spiral arm in IC 342 has an exponential nature; the spi- 
ral arm extends almost linearly on the log R-9 diagram 
(jSato 1 120061: iHirota et al.l 120101) . Lines of the constant 
spiral phase {9 = 0°, —13°) are indicated in the figure for 
comparison. The comparison of the distribution of the 
GMCs with respect to the spiral phases shows that the 
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Figure 6. (a) AV — R plot for the identified GMCs. Crosses in- 
dicate the values corrected with the Gaussian correction, and gray 
filled circles indicate t he uncorrected values . Broken line indicates 
the scaling relation of lSolomon et aLl 1 119871) . For comparison, the 
scaling relations with different coefficients (multiplied by 1.5 and 
1/1.5, respectively) are indicated with dash-dotted lines, (b) Same 
as (a) but for Mco ~ (c) Same as (a) but for Mqq — AV. (d) 
Same as (a) but for Mvir — Mco- 



"wHII" GMCs are located downstream of the "woHII" 
GMCs. If the line oi 6 — —13° is taken as the parti- 
tion line, the GMCs upstream and downstream of the 
line show the difference in terms of the star formation 
activity. The difference in the star formation activity 
within each GMC seems to be well related with the po- 
sition with respect to the spiral arm; while all of the 
clouds upstream of the line are "woHII" GMCs, most 
of the clouds downstream of the line are "wHII" GMCs. 
This is in accordance with the prediction of density wave 
induced/regulated star formation. Hereafter, we will in- 
spect the change of the molecular gas properties accord- 
ing to the classification defined. 



3.6. Line ratio 

Because of the low critical densities of the ^^CO 
(1-0) and the ^^CQ (1-0) lines (nn^ < 10^ cm'^) and 
the large difference of the optical depth between both 
the lines, the line ratio between both the lines (/(^'^CO 
(l~0))//(i2CO (1-0)) = i?i3/i2) could be used as a 
probe which distinguishes discrete cloud structures such 
as ridge of GMCs from diffuse cloud envelopes. It is 
widely known from the observations of Galactic molec- 
ular clouds that while j^i3/i2 is high at the center of 
discrete GMCs (1/3.7; IPolk et al.l fl^ it is low (1/10- 
1/20) at the peripheral regions of GMCs (jSakamoto et al.l 
Il994f) and in small translucent clouds or diffuse high- 
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Table 2 

Basic Molecular Cloud Properties 



ID 


Position'' 


^LSR 


peak 


R 


Ay 








(arcsec, 


arcsec) 


(km s-l) 


(Jy beam~^) 


(pc) 


(km s^-"-) 


(106 Mo) 


(106 Mq) 


1 


16.8 


-17.7 


87.2 


0.47 


62.0 


12.4 


3.58 


1.81 


2 


25.0 


-23.8 


89.5 


0.47 


57.2 


11.6 


3.32 


1.47 


3 


15.7 


1.8 


92.6 


0.40 


44.3 


8.9 


1.32 


0.66 


4 


12.6 


15.9 


93.1 


0.44 


54.5 


9.2 


1.84 


0.87 


5 


-5.3 


17.0 


100.5 


0.41 


62.0 


14.1 


3.66 


2.34 


6 


-0.8 


11.0 


97.0 


0.32 


32.4 


14.6 


1.37 


1.31 


7 


23.1 


-51.8 


79.4 


0.62 


37.2 


9.4 


2.42 


0.63 


8 


12.1 


-7.6 


90.9 


0.28 


69.5 


12.1 


2.39 


1.91 


9 


0.4 


4.8 


87.1 


0.27 


51.5 


13.0 


1.46 


1.66 


10 


7.7 


-0.9 


87.2 


0.27 


45.6 


10.9 


1.46 


1.03 


11 


-4.1 


30.9 


96.3 


0.60 


66.5 


10.2 


3.98 


1.31 


12 


2.7 


15.3 


95.4 


0.29 


42.6 


15.2 


1.80 


1.87 


13 


10.3 


-17.1 


96.6 


0.24 


52.0 


16.0 


1.82 


2.52 


14 


23.0 


-38.0 


84.1 


0.29 


38.9 


9.8 


0.83 


0.71 


15 


0.9 


-14.4 


92.1 


0.24 


65.3 


11.1 


1.34 


1.53 


16 


23.6 


-33.1 


87.9 


0.26 


43.5 


14.8 


0.87 


1.80 


17 


1.4 


20.9 


94.3 


0.31 


55.9 


12.8 


2.11 


1.74 


18 


-7.3 


8.9 


97.3 


0.24 


44.9 


13.4 


1.04 


1.52 


19 


13.7 


-22.4 


93.2 


0.22 


39.9 


16.2 


0.99 


1.98 


20 


6.4 


8.1 


94.1 


0.20 


52.7 


12.3 


1.10 


1.50 


21 


1.3 


-21.8 


93.4 


0.21 


35.1 


9.1 


0.57 


0.55 


22 


-9.7 


-27.1 


90.8 


0.29 


40.7 


10.8 


0.92 


0.89 


23 


-20.3 


5.0 


77.2 


0.26 


41.5 


13.7 


1.02 


1.48 


24 


-14.1 


9.7 


86.2 


0.23 


42.2 


9.8 


0.94 


0.77 


25 


-5.2 


-4.2 


88.5 


0.18 


62.2 


14.9 


1.33 


2.62 


26 


31.4 


-8.4 


92.9 


0.18 


38.1 


9.6 


0.38 


0.67 
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Figure 7. Histogram of the S/im flux measured within each cloud 
boundary of the " wHH" GMCs (gray area) and the " woHH" GMCs 
(white area). 



latitude clo uds fBlitz fc Starl3ll986tlMagnani et al.|[l985l: 
|KnaB2_&iiowcrs 1988). 

Previous extended mapping of IC 342 in the ^'^CO (1-0) 
found that there exist variations of i?i3/i2 in the disk 

^13/12 ■was found to be 



of IC 342 (IHirota et al.llMol) 

low (~ 0.1) at the center of the GMA and both the bar 
ends compared to the other disk regions (0.14-0.2). It 
was also found that the star formation activity in such 
low i?i3/i2 regions are low compared to the high R13/12 
regions. The spatial relation between R13/12 and the star 
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Figure 8. Locations of the identified GMCs superposed on con- 
tour map of the velocity integrated ^^CO (1-0) image. Contour 
levels of the ^^CO images are same as in Figure [2] Crosses and 
circles indicate the locations and the sizes of the GMCs. Radii 
of the circles are taken as an effective radius of the GMCs. Red 
markers represent the GMCs associated with H II region ("wHII" 
clouds), while blue markers represent the GMCs without H II re- 
gion ("woHII" clouds). Lines of constant spiral phase are indicated 
with a dashed line (ip = 0°) and a solid line (</> = —13°). 



formation su ggests that unlike in the starburst galactic 
centers (e.g.. lPaglione et aI]|200lD . the temperature vari- 
ation is not a main cause of the i?i3/i2 variations in the 
disk region. Likely explanation of the i?i3/i2 variations 



in the disk is that R 



13/12 



reflects the fraction of the dif- 
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fuse molecular component within the beam. The diffuse 
molecular component refers to molecular cloud with low 
column and volume densities. 

A map of -R13/12 was made with the following procedures. 
First, to improve the S/N, both the ^^CO and ^^CO data 
cubes were smoothed to 5" resolution 80 pc). Next, 
masking data were made from the ^^CO data cube fol- 
lowing the same procedure described in Section 13.21 and 
were applied to the ^■^CO data cubes. Validity of the 
mask was checked with the same procedure performed in 
Section [3?2] Finally, the R13/12 value for each point was 
calculated using the masked data. 

Figure |9] shows the map of R13/12 compared with the 
^^CO (1-0) map. -R13/12 varies over the observed fields 
from ~ 0.06 to ^ 0.25. An error of the ratio was esti- 
mated following the propagation of the error and was ~ 
0.01 at the brightest ^^CO peaks (Point 1-3) and typ- 
ically below 0.02 at the rest. -R13/12 is rather low on 
the northern side compared to the southern side. Es- 
pecially, at the center of the GMA, R13/12 marks low 
value (~ 0.06 ± 0.01), whi ch is by far small compared to 
the ridge of GMCs a/3.7. [PoIketaI1[T988l) and cl ose to 
the values found in diffuse clouds ~(Blitz & StarSlilll; 
[Magnani et al.l 11985 : Knapp & Bowers 1988). On the 
other hand, southern clouds show higher R13/12 (^ 0.16 
± 0.01), suggesting discrete nature of the clouds. 
Comparing the distribution of the "wHII" and "woHII" 
clouds with the observed i?i3/i2 variation shows that 
there is a possible tendency that while -R13/12 is lower 
in the "woHII" GMCs, it is higher in the "wHII" GMCs. 
The configurations of the R13/12 and star formation dis- 
tributions exclude temperature variation as a primary 
factor for the R13/12 variations. The configurations also 
likely exclude abundance variation induced by selective 
photo dissociation of rarer isotopic species. As the ob- 
served size scale of the R13/12 variation is an order of 
lOOpc, enrichment of ^^C through stellar processing con- 
flicts in terms of timescale. Thus, the variation of the 
^•^CO/^^C O abundance ratio is also unlikely. As sug- 
gested by iHirota et al.l (|2010[ ) , the most likely explana- 
tion for the i?i3/i2 variation in the disk of IC 342 is the 
variation of fraction of the diffuse molecular component. 
Formation of dense cloud core (tt-Hs ~ 10^ cm~^) is re- 
quired for the onset of star formation. Very low -R13/12 
values found at the locations of some of the "woHII" 
GMCs (i?i3/i2 < 0.1) suggest that those "woHII" GMCs 
contain large fraction of diffuse cloud components and 
thus lower fraction of cloud cores compared to the 
"wHII" GMCs. 

4. DISCUSSION 

4.1. Change of Cloud Property Across the Spiral Arm 

The offset between the molecular gas traced by 
the ^^CO(l-O) emission and massive star forming re- 
gions is often seen in grand-design spiral gala xies (e.g., 
IVogel et al.l [T988: Rand et"an[T99a lRa^[l995h . Typical 
extent of the spatial offsets between the two components 
are typically 100-300 pc. However, most of the previous 
observations for such grand-design spiral galaxies, which 
are beyond the Local Group of galaxies, were at resolu- 
tions of 200-1000 pc and were insufficient to investigate 
the variation of the molecular gas properties within these 



0.00 0.05 gjfo'izco 




3h47m15s 10s 5s 47m0s 

Right Ascension (J2000) 



Figure 9. Color representation of /?i3/i2 superposed by contour 
map of ^^CO (1-0) image. Contour levels are 1, 3, 5, and 7 times 
1.5 Jy beam"^ km s~^. Both images were made from the data 
cube smoothed to 5" resolution. Locations of the "wHII" and 
the "woHII" GMCs are indicated with filled circles and crosses, 
respectively. Errors of the R13/12 are typically ~ 0.01 at the bright 
peaks and ~ 0.02 at the rest regions. 

offsets. The NMA-|-45m combined data provide a unique 
opportunity to look into the process involved there with 
a resolution of 50pc and an improved sensitivity to the 
extended diffuse emission. 

Figure [10] shows scatter plots of the properties of clouds 
for each category. The R~ AV plot (Figure rTOT a)) shows 
that "woHII" GMCs have larger line width compared to 
the "wHII" GMCs with the same size. The Kolmogorov- 
Smirnov test was performed to the ratio between R and 
AV with a null hypothesis of both categories being ex- 
tracted from the same parent sample. The null hypothe- 
sis was rejected with a p- value of ^ 0.004. A similar trend 
is also seen in the AV-Mco plot fFigurelTUbV Again, the 
Kolmogorov-Smirnov test was performed and indicated 
that there is a significant difference between the both cat- 
egories with a p-value of ~ 0.001. Another difference of 
properties between the "wHII" and the "woHII" GMCs 
is seen in the Mco-M^n- plot (Figure fTUTc')). There is 
a tendency that the "woHII" GMCs have a larger virial 
mass to luminosity mass ratio than the " wHII" GMCs {p 
^ 0.001), suggesting that "woHII" clouds are less gravi- 
tationally bound compared to the "wHII" clouds. 
Figure [Tl] shows histograms of the basic properties of 
the clouds (radius, line width, and mass). While a little 
difference between the two categories is seen in the ra- 
dius histogram, apparent differences are seen in the mass 
and the line width histograms. As we have seen in the 
previous scatter plots, the line width distribution of the 
"woHII" GMCs seems to be different from the "wHII" 
GMCs (p ~ 0.017, for the Kolmogorov-Smirnov test). 
The median line width of the "woHII" GMCs is - 13.4 
km s~^ and is larger than that of the "wHII" GMCs 
10.2 km s~^). Two cloud categories also show differ- 
ence in the mass histogram {p ~ 0.011). The mass of 
the "woHII" GMCs is concentrated around the median 
value of 1.3 x lO^M©, except for the cloud- 5 with the 
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mass of ~ 3.7 xlO^ M0, which is located at the center 
of the GMA. On the contrary, the mass distribution of 
the "wHII" GMCs is more widely spread and stretched 
to a higher mass regime, up to 4.0 xlO^M©. Except 
for a cloud-26, aU masses of "wHII" GMCs are larger 
than IO^Mq. The cloud-26, which is located ~ 300 pc 
downstream from the molecular ridge, has mass of 0.4 
X 1O^M0 and is likely being disrupted by feedback from 
associating H II regions. 

The comparison between the cloud properties between 
the two cloud groups (the "wHII" and the "woHII" 
GMCs) indicated significant variation of cloud properties 
according to their associated star forming activity. The 
"woHII" GMCs have larger line width and have smaller 
mass and moreover, are less gravitationally bound com- 
pared to the "wHII" GMCs. Spatial distribution of the 
clouds indicate that the "wHII" GMCs are located down- 
stream of the "woHII" GMCs, on average. These facts 
suggest that cloud properties do change by crossing the 
spiral arm. 

4.2. Effect of Gas Velocity Dispersion 

The comparison between the virial mass and the 
luminosity-based mass made in the previous subsection 
indicated that the "woHII" GMCs are loosely bound 
compared to the "wHII". The degree of gravitational 
binding can be expresse d with a virial parameter, (ofvii- 
iBertoldi fc McKeg|1992[ ) , which is expressed by 

ttvir = Afvir/Mco- (6) 

As ttvir decreases, cloud gets more unstable against grav- 
itational collapse. In this subsection, we will utilize ckvir 
to investigate which of the cloud properties is responsible 
for the observed variation of the gravitational instability 
of the clouds. 

Figure [T^ shows the scatter plots between the basic cloud 
properties and avir- To see whether there exist a sig- 
nificant relation between the virial parameter avir with 
other cloud properties, Kendall's r and associated prob- 
ability of chance correlation were calculated. Significant 
correlation was only detected between the virial param- 
eter ttvir and the line width (p ^ 6.0 x 10~^). GMCs in 
the Galactic center exhibit similar correlation between 
the line width a nd the degree of gravitational binding 
(|Oka et al.ll2001l ). The correlation suggests that the de- 
gree of gravitational binding of the clouds is mainly de- 
termined by their internal turbulent motion. 
The "wHII" GMCs have smaller velocity dispersion and 
avir, and located downstream compared to the "woHII" 
GMCs. Similar relation between the star formati on ac- 
tivity and the line width was pointed out by Kohn o~et al.l 
([1999) . They found an anticorrelation between the dense 
molecular gas fraction and the line width in the center 
of the barred galaxy NGC 6951, with a spatial resolu- 
tion of ~ 400pc. The results presented here suggest that 
the anticorrelation between the star formation and the 
line width holds down to the size scale of CMC. It is 
tempting to consider that dissipation of turbulent mo- 
tion inside the GMCs occurs by crossing the spiral arm, 
and leads to the sub sequent massive star formation. 
We must note that iScoville et al.l (|1986[ ) indicated the 
relation between the star formation activity and the line 
width, but in an opposite sense; GMCs smaller than 40pc 



show that the GMCs associated with H II regions have 
larger line width compared to the GMCs without H II 
regions. The difference in the results is likely due to 
the difference in the observed siz e scales. Since the sp a- 
tial resolution of the data used bv 'Sc oville et al.l ()1986( ) is 
finer than this work (down to pc scale), the effects of local 
phenomena, such as shock induced by cloud-cloud col- 
lision and radiation feed backs from H II regions, might 
be dominant in the previous result. While, on the other 
hand, the resolution of this study is comparable to the 
size of the GMCs and the derived properties should be 
reflecting the 'global' properties of clouds rather than the 
internal phenomena within each cloud. 

4.3. Evolution of GMC in the Spiral Arm 

The results presented in the previous subsections il- 
lustrated the variations of the cloud properties across the 
spiral arm in IC 342. The "wHII" (downstream) GMCs 
are more virialized and massive compared to the " woHII" 
(upstream) GMCs. The cloud groups also differ from 
each other in that the velocity dispersion of the "wHII" 
GMCs is larger than that of the "woHII" GMCs. In this 
subsection, we will discuss the implications of such vari- 
ations of the cloud properties. 

First question is from what kind of material the clouds 
identified in this study were made? There are two classes 
of concepts about the formation of GMCs in the spiral 
arm. One is that GMCs are formed from H I clouds in the 
inter-arm and another is that GMCs are built-up from 
smaller clouds in the inter-arm through agglomeration. 
As we have seen before in Suction [3321 the surface mass 
density of H I is only ~ 2Mq pc^^ within the observed 
region. On the other hand, the surface mass density of 
the molecular gas is estimated to be ^ 3OAf0 pc~^ from 
the single-dish ^^CO data. So, it is natural to consider 
that the GMCs in the spiral arm are formed from the 
pre-existing molecular clouds. 

As the extent of the observed region is limited, we are un- 
able to tell the nature of the inter-arm molecular clouds 
from the observed data alone. Instead, we refer to the 
observations of molecular clouds in the Milky Way. It 
is known that GMC with mass ab ove 10^ is rarel y 
seen in the i nter-a rm region (e.g., iSanders et all 119851 ). 
iHever et al.l ([1998D indicated that inter-arm molecular 
clouds are smaller than 3 x IO'^Mq and not virialized 
in the outer Galaxy. Even in the molecular-rich in- 
ner Galaxy, there are substantial contributions of dif- 
fuse mole cular emiss i on, suggesting the existe nce of small 
clouds. ([Polk et al.lll988HChiar et al.lll994D . It seems 
natural to assume that inter-arm molecular clouds are 
dominated by the small molecular clouds. 
Second question is what kind of mechanism could 
explain the observed variations the of cloud proper- 
ties? If the inter-arm clouds in IC 342 are small, dif- 
fuse clouds, the mechanism for building up the GMCs 
as massive as IO^Mq must exist. The timescale 
for building up should be comparable to the arm 
crossing tirne in this region (tarm ^ 3 x lO^yr, 
iHirota et al.ll2010l ). The random coalescence of molecu- 
lar clouds is proposed as a forma tion mec hanism of GMC 
(iKwan fc V aldcs 1983; Casoli & Combesl[l 982: To misaka 
19841 [Roberts fc Stewartiil987i: iElmegreenlil990i: iDobbs 
2008[ ). Early calculations indicated that it may take 
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Figure 10. a) AV — i? plot of the identified GMCs. Red and blue crosses indicate the "wHII" GMCs and the "woHII" GMCs, respectively. 
Broken line indicates the scaling relation of Solomon et al. (1987). For comparison, the scaling relations with different coefficients (multiplied 
by 1.5 and 1/1.5, respectively) are indicated with dash-dot lines, b) same as a), but for Mqq — t^V . c) same as a), but for -A/vir — Afco- 



about 10* yr, which is far longer than <arm- However, 
if the spiral density waves are included, it turns out that 
the formation of the GMC as massive as several times 
IO^A/q within a time comparable with iarm is possible. 
Large-scale instability induced by the self-gravity of gas 
shall further make the agglomeration p rocess efhcient 
(e.g., IBalbus fc Co^fToMlDobbsl [20081 ). 
If the random coalescence model is adopted, the observed 
nature of the "wHII" and the "woHII" GMCs could 
be understood naturally. GMCs were made from the 
small molecular cloud in the inter-arm through agglom- 
eration. The end product of the cloud growth will be the 
"wHir'C'post arm") GMCs, which have masses of (1-3) 
X IO^Mq (an exception is the cloud-26). As the masses of 
the "woHII" GMCs are around lO^Af0 and smaller com- 
pared to the " wHII" GMCs, it may be natural to consider 
that the " woHII" GMCs are at the intermediate stage of 
the cloud growth by coalescence. The relatively low peak 
temperatures of the "woHII" GMCs (Tmb = 2-3K) and 



the low degree of gravitational binding (avir = 1~3) are 
in agreement with the notion that the "woHII" GMCs 
are collective of smaller clouds, which are yet not virial- 
ized. 

Above interpretation may be in agreement with the re- 
cent results of observations of the grand-design spiral 
galaxy M51. iKoda et "all (|2009| ) performed full aperture 
observation of M51 in -'^^CO (1-0) with a spatial res- 
olution of ^ 200 pc and indicated that while massive 
GMAs are preferentially found in the spiral arm, GMCs 
with mass below 10^ are uniforrn l y dist ributed over 
the molecular gas disk. lEgusa et al.l (I20TTI ) carried out 
CO observations toward the selected region of M51 with 
a spatial resolution of 30pc, comparable with ours, but 
not corrected for missing-flux, and found that clumps 
with mass of 10^-10® Mq are likely preferentially lo- 
cated on the downstream side of the spiral arm. The 
combined data presented here revealed while rather dif- 
fuse " woHII" clouds are the present at the upstream side, 
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Figure 11. Histogram of the basic properties of the GMCs, 
namely, deconvolved effective radius, deconvolved line width, and 
luminosity mass. Histograms for the "wHH" GMCs, the "woHII" 
GMCs, and all GMCs are indicated with red solid, blue solid, and 
dotted lines, respectively. Vertical dashed lines in red and blue 
indicate the median value of each category. 
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discrete clouds are seen on the downstream side. While 
the buildup of GMAs by coagulati on of GMCs in th e 
inter-arm of M51 was suggested by iKoda et al.l ()2009D . 
the low i?i3/i2 in the center of the GMA and presence of 
diffuse clouds in the upstream of the spiral arm implies 
coalescence of diffuse clouds are much important in the 



arm of IC 342. 

Another important aspect of the variation of cloud prop- 
erties is the decrease of line width. Since molecular 
clouds collide inelastically, part of the kinetic energy 
should be transf ormed into thermal energy and radiated. 
iTomisakal (|T987^ carried out the numerical simulation of 
collisional clouds incorporating formation and destruc- 
tion of GMCs, the loss of random velocity due to inelastic 
cloud collision, and the energy input by the star forma- 
tion within GMCs. The results were that after the com- 
pression of cloud by the shock, random velocity of clouds 
decreases and density increases. These are in agreement 
with the results presented in previous subsections. 
The scenario of cloud growth would be summarized as 
follows. First, the inter-arm molecular clouds are col- 
lected because of the convergence of stream line before 
the spiral arm. Due to this convergence of stream line 
and self-gravity, clouds will stick and form the cloud col- 
lective with masses of lO^M©, although they are not fully 
virialized yet and loosely bound. These cloud collectives 
are observed as the "woHII" GMCs. Once massive CMC 
is formed, its large gravitational cross section accelerates 
the coalescence and swallows diffuse clouds at the enve- 
lope. The excess kinetic energy will be dissipated due to 
inelastic collision and virialized GMCs with mass larger 
than IO^Mq are formed (the "wHII" GMCs"). 

4.4. Condition of Massive Star Formation 

The closeness between the star forming regions and 
the "wHII" GMCs suggests that these clouds are real 
progenitors of the associated H II regions. This can also 
be justified by the consideration of the age of H II regions. 
Sizes of the H II regions seen in the Ha image have at 
most a diameter of ~ 100 pc. Assuming sound speed in 
the ionized gas as c ~ 10 km s~^, an upper limit on the 
age of the H II regions could be estimated as D/2c = 5 
Myr. As the age of H II regions estimated is much 
shorter than the arm crossing time (tarm 3 x IO^Mq, 
), the H II regions must have been born within the as- 
sociated "wHII" GMCs on site. The relatively young 
age of the H II regions compared to iarm gives us an- 
other important implication; if the "woHII" GMCs are 
able to initiate massive star formation, we should see the 
symptoms of massive star formation inside the " woHII" 
GMCs. Nevertheless, little 8/im emission could be seen 
in the "woHII" GMCs. This implies that there must be 
a condition for massive star formation, which is fulfilled 
by the "wHII" GMCs but not by the "woHII" GMCs. 
As we have seen in Section 14. 3[ collisional coagulation 
of clouds seems to be responsible for the evolution of 
the GMCs. One of the possible mechanism for ini- 
tiating massive star formation is cloud-cloud collision 
(|Scoville et al.lll986l : [^i2000( ). In this picture, molec- 
ular clouds are compressed by the shock at the col- 
liding surface and sub sequently dense cores are formed 
(jKimura fc Tosal [19960 . However, at the center of the 
GMA where large velocity gradient due to streaming mo- 
tions is observed and enhanced rate of collision is ex- 
pected, little or no sign of star formation is seen. This 
implies that for the onset of massive star formation, not 
only cloud collision but also some other conditions may 
be required. Another interpretation as following might 
arise. Namely, after star formation is initiated inside the 
"woHH" clouds, star formation proceeds rather slowly 
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and become only visible at the "wHII" clouds. How- 
ever, this requires slow progress of star formation, with 
a timescale comparable to the arm-crossing time men- 
tioned above. Even if this is the case, the progress of 
the star formation shall be regulated by the evolution 
of GMCs because of the far longer timescale of arm- 
crossing time compared to dynamical time scale in each 
GMC (~ 4x lOV, e.g., iPringle et aLll200H) . Both inter- 
pretations requires the evolution of the cloud properties 
which enable triggering of star formation to work effec- 
tively. 

Figure [T51 shows a scatter plot of the area-averaged 8/im 
luminosity measured within each cloud boundary and the 
virial parameter Ovir for each cloud. There is a tendency 
that 8/im luminosity, which is proportional to the star 
formation rate, is higher in the range of avir < 1 com- 
pared to in the range of avir > 1- It was shown in Section 
14.21 that the degree of binding of the GMCs is mainly 
determined by its line width. These facts suggest that 
virialization of GMC achieved by the dissipation of tur- 
bulent motion is one of the necessary conditions for the 
onset of massive star formation. 

Recent studies emphasize the importance of turbu- 
lence in both regulating and the triggering star forma- 
tion. In the theory of turbulence regulated star for- 
mation, although dense cores are formed by the shock 
driven by supersonic turbulence inside a GMC (e.g . , 
IPadoan fc Nordlundl [200l IMac Low fc KlessenI [200l . 
limited part of the GMC can be dense enough to form 
cloud cores since the greater part of the GMC are sup - 
ported by the exces s kine tic energy ()Elmegreenl I2Q02D . 
iKrumholz fc McKed ()2005[ ) presented semianalytical cal- 
culation using the theory and showed that amount of the 
excess kinetic energy inside the GMC, which is propor- 
tional to avir, controls the production rate of cloud cores, 
and thus star formation rate. This point is in agreement 
with our results that dissipation of turbulence is required 
for the onset of massive star formation. After GMCs 
have dissipated excess kinetic energy, it can initiate star 
formation, whatever the trigger may be. 

5. SUMMARY 

Results of the full aperture observations of the north- 
eastern spiral arm segments of IC 342 in ^^CO (1-0) and 
^■^CO (1-0) made with a spatial resolution of ~ 50 pc 
were presented. Though the spatial resolution of the ob- 
servations is not high enough to resolve internal struc- 
ture of GMCs, it is capable of investigating the global 
property of massive GMCs. The observations cover the 
Ikpc X 1.5kpc region which contain a GMA with mass of 
~ 10^ Mq, where spatial offsets between ^^CO and star- 
forming regions exist. The NMA data were combined 
with the 45m telescope data to image the extended CO 
distribution and investigate whether and how the prop- 
erties of molecular clouds change by crossing the spiral 
arm. The results are summarized as follows. 

1. The spiral arms were resolved into a number of 
clouds which have size, temperature, and surface 
mass density comparable to the massive GMCs in 
the Galaxy. While the ^^CO (1-0) distribution 
is concentrated in a narrow ridge in the southern 
half of the observation field, it is more smoothly 
distributed in the northern half, where the GMA 
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Figure 13. Scatter plot of area-averaged 8/im luminosity mea- 
sured within each cloud boundary and the viral parameter (ctvir) 
for each cloud. Filled circles indicate the "wHII" GMCs, open 
diamonds the "woHII" GMCs, and filled gray circle indicates the 
cloud-26. 



exists. Comparison with the ^'^CO (1-0) image 
showed that at the center of the GMA only lit- 
tle ^'^CO (1-0) emission is detected despite of its 
strong ^^CO emission. 

2. Comparison with star formation tracers indicated 
that while some of the clouds are closely associated 
with star forming regions, the rest of the clouds 
show little or no sign of massive star formation. 
The identified clouds were divided into two cat- 
egories according to whether they are associated 
with star formation activity or not. 

3. Twenty-sixclouds with masses of (0.4-3.7) x 10^ 
Mq were identified from the ^^CO (1-0) data cube. 
The size and the line width of the clouds were com- 
parable to GMCs in the Milky Way. The iden- 
tified clouds also followed the linewidth-size rela- 
tion and the mass-lincwidth relation of the Galac- 
tic GMCs, suggesting that the identified clouds are 
resemblance of Galactic GMCs. 

4. The identified clouds (GMCs) were divided into 
two categories according to whether they are asso- 
ciated with star-forming regions or not. Compar- 
ison between both categories showed that clouds 
which are associated with star-forming regions 
("wHII" GMCs) are more virialized and massive 
compared to the clouds which show little or no sign 
of star formation ("woHII" GMCs). Moreover, the 
"woHII" GMCs have larger line width compared 
to the "wHII GMCs". As the "wHII" GMCs are 
located downstream of the "woHII" GMCs, it is 
concluded that properties of the GMCs do change 
by crossing the spiral arm. 

5. The line ratio, R13/12, varies from 0.06 ± 0.01 in 
the center of the GMA where little or no sign of 
star formation is seen to ~ 0.16 in the discrete 
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GMCs which arc associated with H II regions. On 
the whole, the "woHII" GMCs likely show a low 
-R13/12 value which suggests a larger fraction of dif- 
fuse molecular cloud components than the "wHII" 
GMCs. 

6. To see what parameter of the GMC controls 
degree of gravitational boundness of the GMC, 
correlation plots between the virial parameter, 
ttvir = Myir/Mco, and other basic GMC prop- 
erties, namely, radius, line width, and mass were 
made. The linewidth-avir plot indicated most 
strong correlation implying that the dissipation of 
turbulent motion controls the boundness of the 
GMCs. 

7. Random coagulation of the pre-existing small 
clouds can explain the growth of the mass and the 
timescale involved in the change. CoUisional coa- 
lescence involved in random coagulation may also 
explain the observed decrease of line width, as the 
kinetic energy of the clouds will be dissipated by 
inelastic cloud collision. 

8. The closeness of H II regions and the " wHII" GMCs 
imply that those GMCs are real progenitors of the 
associated H II regions. The rough estimate of the 
age of H II regions gives ^ 5 Myr, which is well 
below the arm-crossing time. This implies if the 
"woHII" GMCs are able to initiate star formation, 
the symptoms should have seen. So, the sparseness 
of star formation in the "woHII" GMCs suggests 
they do not meet the condition for massive star 
formation. 

Plot of area-averaged 8/xm luminosity and virial pa- 
rameter (avir) shows a tendency that the 8/xm lumi- 
nosity, which is proportional to the star formation 
rate, rises in the range of a^ir < 1- This might im- 
ply that dissipation of turbulence controls not only 
the evolution of GMCs, but also the probability of 
massive star formation. 
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